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Abstract 

We examine e + e~ — » D~ D* + and D* s ~ Df interactions at 4170 MeV using the CLEO-c detector 
in order to measure the decay constant f n + with good precision. Previously our measurements were 

U s 

substantially higher than the most precise lattice based QCD calculation of (241±3) MeV. Here 
we use the Df — > i + u channel, where the i + designates either a fi + or a r + , when the r + — > 7r + F. 
Analyzing both modes independently, we determine B(D+ — > ^ + z>) = (0.565 ±0.045 ±0.017)%, and 
B(Df — > r+z^) = (6.42 ± 0.81 ± 0.18)%. We also analyze them simultaneously to find an effective 
value of 6 cff (L>+ -► /x+i/) = (0.591±0.037±0.018)% and / D + = (263.3±8.2±3.9) MeV. Combining 
with the CLEO-c value determined independently using Df — > r + zv, r + — > e + i^ decays, we extract 
/ D + = (259.5 ± 6.6 ± 3.1) MeV. Combining with our previous determination of B{D + — > /U + ^), we 
extract the ratio f D +/fo+ = 1.26 ± 0.06 ± 0.02. No evidence is found for a CP asymmetry between 
r(D+ — > /i + ^) and r(D~ — > ^~z^); specifically the fractional difference in rates is measured to be 
(4.8±6.1)%. Finally we find B(D+ -► e+zv) < 1.2 x 10" 4 at 90% confidence level. 
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I. INTRODUCTION 



We discuss here an improved measurement of the width of the purely leptonic decay 
— > £ + u, when the £ + is either a /x + or a r + , when the latter decays into a tt + V [1]. In a 
companion article we report an improved measurement of the decay width for — > r + z/, 
where r + — > e + vv. 

In the Standard Model (SM) these decays are described by the annihilation of the initial 
quark- ant iquark pair into a virtual W + that materializes as a l^v pair; the process is shown 
in Fig. [TJ The decay rate is given by [3| 



r»2 



where M D + is the Df mass, rri£ is the mass of the charged final state lepton, Gp is the 
Fermi coupling constant, and \V CS \ is a Cabibbo-Kobayashi-Maskawa matrix element with a 
value we take equal to \V u d\ of 0.97418(26) pj, and / D + is the "decay constant," a parameter 
related to the overlap of the heavy and light quark wave-functions at zero spatial separation. 




FIG. 1: The decay diagram for — > 

The SM decay rate then is predicted using a theoretical calculation of the decay constant. 
Two calculations have been carried out using unquenched lattice quantum-chromodynamics 
(LQCD). Aubin et al. find f D + = (249 ± 3 ± 16) MeV while a more recent calculation 
of Follana et al. gives (241 ± 3) MeV (a|. The latter calculation is more than three standard 
deviations lower than the average of previous CLEO and Belle measurements [3]. 

Dobrescu and Kronfeld have proposed three models based on physics beyond the SM that 
are consistent with known data and could possibly explain the difference. One is a charged 
Higgs model and the other two involve different manifestations of leptoquarks The recent 
CLEO measurement of fn+ = (205.8 ± 8.5 ± 2.5) MeV, is consistent with both the Aubin et 
al. and Follana et al. predictions, of (201 ±3 ± 17) MeV and (208 ±4) MeV, respectively Q. 

It is particularly important to understand if the discrepancy in the D s case is due to 
physics beyond the SM, a faulty theoretical calculation, or to an unlikely measurement fluc- 
tuation. We note that precise information on the size of Cabibbo-Kobayashi-Maskawa matrix 
elements extracted from B — B mixing measurements requires theoretical input on the "de- 
cay constants" for B d and B s mesons or their ratio, JbJ f'B d [10]. Although the calculations 
in the B and D systems are not the same, many of the techniques used are common, and a 
discrepancy in the charm system at a minimum, does not give confidence in the theoretical 
predictions for the B system. In this paper we present an updated measurement of f D + with 
much improved precision. 
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Akeroyd predicts that the presence of a charged Higgs boson would suppress fo 
There is however the possibility, not considered by Akeroyd, that it is the charm quark that 
is responsible for a NP contribution not the s quark [8j. In that case the relative change 
would be similar in D + and D£ decays. 

We can also measure the ratio of decay rates to different leptons, and the SM predictions 
then are fixed only by well-known masses. For example, for t + u to 
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Using measured masses [121 ]. this expression yields a value of 9.76 with a small error. Any 
deviation in R from the value predicted by Eq. (|2D would be a manifestation of physics beyond 
the SM. This could occur if any other charged intermediate boson existed that affected the 
decay rate differently than mass-squared. Then the couplings would be different for muons 
and r's. This would be a clear violation of lepton universality (l3| . 

Most other measurements of f D + have been hampered by a lack of statistical precision, 
and relatively large systematic errors Il5l . Il6l . 1 1 71 . Il8l . [19]. One large systematic error 



source has been the lack of knowledge of the absolute branching fraction of the normalization 
channel, usually — > <f)Tr + 2(|. The results we report here will not have this limitation, 



nor did our previous measurement [l|, nor did the Belle measurement [211 ] . 

In both /i + z/ and r + u D s decays the charged lepton must be produced with the wrong 
helicity because the D s is a spin-0 particle, and the final state consists of a naturally left- 
handed spin-1/2 neutrino and a naturally right-handed spin-1/2 anti-lepton. Because the 
r + has a mass close to that of the Df, the helicity suppression is broken with respect to the 
/x + decay, but there is an additional large phase space suppression. Because of the helicity 
suppression in fi + u the radiative process 7/i + z/ may have a significant rate. Dobrescu and 
Kronfeld, however, estimate this process is only 1% of the lowest order mechanism, for 
photon momenta below 300 MeV, which is relevant range for this analysis. We include this 
radiative correction in what follows 

(There is no correction for the t + u final state.) 
The data will also be corrected for final state radiation of the muon, as our Monte Carlo 
simulation incorporates this effect [23| . 



II. EXPERIMENTAL METHOD 
A. Selection of D s Candidates 

The CLEO-c detector 24j is equipped to measure the momenta and directions of charged 
particles, identify them using specific ionization (dE/dx) and Cherenkov light (RICH) [251 ]. 
detect photons and determine their directions and energies. 

In this study we use 600 pb _1 of data produced in e + e~ collisions using the Cornell Elec- 
tron Storage Ring (CESR) and recorded near a center-of-mass energy (Ecu) of 4.170 GeV. 
At this energy the e + e~ annihilation cross-section into D~D* + + D*~D^" is approximately 
1 nb, while the cross-section for D^D~ is about a factor of 20 smaller. In addition, D mesons 
are produced mostly as D*D*, with a cross-section of ~5 nb, and also in D*D + DD* final 
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states with a cross-section of ~2 nb. The DD cross-section is a relatively small ~0.2 nb 
26[. There also appears to be DD*n production. The underlying light quark "continuum" 
background is about 12 nb. The relatively large cross-sections, relatively large branching 
fractions, and sufficient luminosities allow us to fully reconstruct one D s as a "tag," and 
examine the properties of the other. In this paper we designate the tag as a D~ and ex- 
amine the leptonic decays of the -D+, though in reality we use both charges for tags and 
signals. Track requirements, particle identification, 7r°, rj, and K s selection criteria are the 
same as those described in Ref. j27[, except that we now require a minimum momentum of 
700 MeV/c for a track to be identified using the RICH. 

We also use several resonances that decay via the strong interaction. Here we select 
intervals in invariant mass within ±10 MeV of the known mass for rf — > 7T + 7T~i], ±20 MeV 
of the known mass for rj' -> p°7, ±10 MeV for <p -> K+K~, ±100 MeV for K*° -> K~ir + , 
and ±150 MeV for p~ — > tt^tt or p° — > 7t + tt~ . 

We reconstruct tags from either directly produced D s mesons or those that result from 
the decay of a D*. The beam constrained mass, mec, is formed by using the beam energy 
to construct the D* candidate mass via the formula 



m BC 



(3) 



where i runs over all the final state particles. If we ignore the photon from the D* — > 7_D S 
decay, and reconstruct the mec distribution, we obtain the distribution from Monte Carlo 
simulation shown in Fig. El The narrow peak occurs when the reconstructed D s does not 
come from the D* decay, but is directly produced. 
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FIG. 2: The beam constrained mass mec from Monte Carlo simulation of e + e~ — ► D~^D* S ~ at an 
Ecm of 4170 MeV. The narrow peak is from the Df and the wider one from D*~ — > jDj . (The 
distributions are not centered at the Df or D* + masses, because the reconstructed particles are 
assumed to have the energy of the beam.) 



Rather than selecting events based on only mec, we first select an interval that accepts 
most of the events, 2.015 < wbc < 2.067 GeV, and examine the invariant mass. Distributions 
from data for the nine tag decay modes we use in this analysis are shown in Fig. [31 Note 
that the resolution in invariant mass is excellent, and the backgrounds not abysmally large, 
at least in these modes. To determine the number of D~ events we fit the invariant mass 
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distributions to the sum of two Gaussians centered at the D s 
as "two-Gaussian." The r.m.s. resolution (a) is defined as 



mass, a function we refer to 



o- = /i<ti + (1 - fi)a 2 , 



(4) 



where a\ and 02 are the individual widths of each of the two Gaussians and fi is the fractional 
area of the first Gaussian. The number of tags in each mode is listed in Table [H There are 
two changes in modes from our previous analysis [1]; instead of <pp~ , we now use K + K~7t~tt°, 
and we have added the 7]'ti~ , 77' — > p°7 mode; the background here is somewhat reduced as 
we apply a cut on the helicity angle, 9 of | cos#| < 0.8, since the p° is polarized as sin 2 6. 
(Here 9 is the angle of the 7r + in the p° rest frame with respect to the p° direction in the 
parent frame.) These changes results in an increase of 20% more tags at the expense of more 
background. 
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FIG. 3: Invariant mass of Dj candidates in the decay modes (a) K + K~tt~ , (b) KsK~, (c) 
r}7r~ (77 — > 77), (d) 7]'ir~ (7/ — > 7r + 7r~7?, 77 — > 77) (e) K + K~7r~7r° (f) 7r + 7r~7r~, (g) K*~~K*°, 
K*~ — > Kg7r~ (K*° — > K + ir~), (h) T/p - , and (i) 77'^^ (77' — > p 7), after requiring the total energy 
of the D~ candidate to be consistent with the beam energy. The curves are fits to two-Gaussian 
signal functions plus a linear background. The dashed vertical lines in (b) indicate the restricted 
fit region. 
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TABLE I: Tagging modes and numbers of signal and background events, within ±17.5 MeV of the 
D~ mass for each mode, determined from two-Gaussian fits to the invariant mass plots, and the 
number of tags in each mode including the 7 from the D* — > 7_D S transition, within an interval 
3.872 < MM* 2 < 4.0 GeV 2 , as determined from fits of the MM* 2 distributions (see text) to a signal 
Crystal Ball function (see text) and two 5th order Chebychev background polynomials. 



Mode 






Invariant Mass 
Signal Background 


MM* 2 
Signal Back 


ground 


K+K-7T- 






26534±274 


25122 16087±373 


39563 


K S K- 






6383±121 


3501 


4215±228 


6297 


r/vr~; 77 -i 


>77 




2993 ± 156 


5050 


2005±145 


5016 


t/tt - ; rf - 


-> 7T + 7r~ 




2293±82 


531 


1647±131 


1565 


K+K~TT' 






11649±754 


78588 


6441±471 


89284 








7374±303 


60321 


5014±402 


43286 


K*-K*°; 


K*~ - 


+ Kfyr-, K*° - 


^ K+TT- 4037±160 


10568 


2352±176 


12088 


rjp-; rj -> 


77, P~ 


- - t^tt 


5700±281 


24444 


3295±425 


24114 


ri'ir-; rf - 


-A, 




3551±202 


19841 


2802±227 


17006 


Sum 






70514 ± 963 


227966 43859±936 


238218 



We list the number of signal events in each mode in Table d by finding the number of 
events within ±17.5 MeV of the D s mass; here we integrate the two-Gaussian PDF over 
the interval. We also include the amount of background in this interval. For ease of further 
analysis we sum all tag modes together, as shown in Fig. HI 

In this analysis we look for two types of events: (1) e + e~ — > D*~Df, and (2) e + e _ — ► 
D~D* + , where our convention is that the tag is denoted by the negative charge state and the 
putative signal fi + u decay by the positive charge state. Thus, we need to detect the photon 
from the D* decay. Therefore, we look for an additional photon candidate in the event that 
satisfies our shower shape requirement. Regardless of whether or not the photon forms a D* 
with the tag, for real D*D S events, the missing mass squared, MM* 2 , recoiling against the 
photon and the D~ tag should peak at the mass-squared. We calculate 

MM* 2 = (E CM - E Ds - E,f - (p CM - p Ds - P 7 ) 2 , (5) 

where Eqm (pcm) is the center-of-mass energy (momentum), Ec a (pdJ is the energy (mo- 
mentum) of the fully reconstructed D~ tag, and E 7 (p 7 ) is the energy (momentum) of the 
additional photon. In performing this calculation we use a kinematic fit that constrains the 
decay products of the D~ to the known D s mass and conserves overall momentum and en- 
ergy. All photon candidates in the event are used, except for those that are decay products 
of the D~ tag candidate. 

The MM* 2 distributions for events in the D~ invariant mass signal region (±17.5 MeV 
from the D s mass) are shown in Fig. [5J In order to find the number of tags used for 
further analysis we preform a two-dimensional binned maximum liklihood fit of the MM* 2 
distribution and the invariant mass distribution in the interval ± 60 MeV from the D s 
mass and 3.5 < MM* 2 < 4.25 GeV 2 . This procedure is improved by having information 
on the shape of the MM* 2 signal function (often called a Probability Distribution Function, 
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1631108-029 




FIG. 4: Invariant mass of Dj candidates summed over all decay modes and fit to a two-Gaussian 
signal shape plus a straight line for the background. The vertical dot-dashed lines indicate the 
±17.5 MeV definition of the signal region. 



or PDF). One possibility is to use the Monte Carlo simulation for this purpose, but that 
would introduce a relatively large systematic uncertainty. Instead, we use our relatively 
large sample of fully reconstructed D S D* events, where we use the same decay modes listed 
in Table [U we find these events and then examine the signal shape in data when one D s is 
ignored. 

To remove background we subtract invariant mass sidebands. Some random photon back- 
ground remains. We remove this by examining the direction of the candidate photon; it 
should be opposite the direction of the D+D~ system for signal. Denning the angle # 7 to 
be that between the three-vector of the candidate photon and the D+D~ system, we require 
that —1.0 < cos(0 7 ) < —0.95. Furthermore we use the events for —0.8 < cos(6 l 7 ) < 1.0 as 
background and subtract them also. The MM* 2 distribution from this sample is shown in 



Fig. [61 The signal is fit to a Crystal Ball function [28|, [29J. The a parameter, that represents 
the width of the distribution, is found to be 0.035±0.001 GeV 2 . We do expect this to vary 
somewhat depending on the final state, but we do not expect the parameters that fix the 
shape of the tail to change, since they depend mostly on initial state radiation, beam energy 
spread, and the properties of photon detection. 

The background has two components, both described by 5th order Chebyshev polyno- 
mials; the first comes from the background under the invariant mass peak, denned by the 
sidebands, and the second is due to multiple photon combinations. In both cases we allow 
the parameters to float. 

We find a total of 43859±936 events within within the interval 3.872 < MM* 2 < 4.0 GeV 2 
and having an invariant mass within ±17.5 MeV of the D s mass, where the total number of 
events is the sum of the yields from the fits to each mode as shown in Table [H An overall 
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MM* 2 (GeV 2 ) 

FIG. 5: The MM* 2 distribution from events with a photon in addition to the D~ tag for the 
modes: (a) K + K~n~, (b) K%K~, (c) t/tt - , (d) rfir~, (e) K+R-tt-tt , (f) vr+vr-vr-, (g) K*~K*°, 
(h) T7/9 - , and (i) ??'vr _ , rj — > 7r + 7r~7 . The curves are fits to Crystal Ball functions and two 5th 
order Chebychev background functions (see text). 



systematic error of 2.0% on the number of tags is assigned by using different functions for 
the description of the backgrounds. If we fix the shape of the multiple photon combinations 
background polynomial to that given by the Monte Carlo, we increase the yield by 1.1%. If 
we use a 4th order polynomial to describe the background of non-D" events, we decrease the 
yield by 1.6%; on the other hand using a 6th order polynomial increases the yield by 1.7%. 
Combining the results of changing both background shapes in quadrature for the worst case, 
we assign a systematic error of ±2.0% on the tag yield. 

We also determine the number of D~ tags when the tag results from the decay D*~ — > 
7-D~ . This procedure reduces the background considerably in determining the tag yield and 
results in counting half of the tags. Adopting this procedure we find a tag yield that is higher 
by (1.2±2.3)% from our nominal procedure, consistent with the assigned systematic error. 
The summed MM* 2 distribution is shown in Fig. [7J 

There is also a small enhancement of 5.2% in our ability to find tags in n + u (or r + z/, 
t + — ► n + v) events (tag bias) as compared with events where the decays generically. We 
determine this correction by using a Monte Carlo simulation for each tag mode independently 
and then average the results based on the known tag fractions. We assign a systematic error 
of 20% giving a correction of (5.2 ± 1.0)%. 
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B. Signal Reconstruction 

We next describe the search for — > fi + u. We select events within the MM* 2 region 
shown in Fig. [7] for further analysis. We note that the limits are rather wide. We use 
this selection because the background in the signal side is rather small and the errors are 
minimized by taking as many tags as possible. 

Candidate events are selected that contain only a single extra track with opposite charge 
to the tag. The track must make an angle >25.8° with respect to the beam line to ensure that 
it is well measured, and in addition we require that there not be any neutral cluster detected 
in the calorimeter, not associated with the tag, with energy greater than 300 MeV (photon 
veto). These cuts are highly effective in reducing backgrounds. The photon energy cut is 
especially useful to reject Df — ► 7r + 7r°, should this mode be significant, and — > ?77r + . 

Since we are searching for events where there is a single missing neutrino, the missing 
mass squared, MM 2 , evaluated by taking into account the observed Dj, and 7 should 
peak at zero; the MM 2 is computed as 

MM 2 = (E CM - E Ds -E 1 - Erf - (p CM - p Da ~ P 7 - P^) 2 , (6) 

where E^ (p M ) are the energy (momentum) of the candidate muon track and all other vari- 
ables are the same as defined in Eq. §5§. 

We also make use of a set of kinematical constraints and fit each event to two hypotheses 
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FIG. 7: The MM* 2 distribution summed over all modes. The curves are fits to the number of 
signal events using the Crystal Ball function and two 5th order Chebychev background functions 
(see text). The vertical lines show the region of events selected for further analysis. 



one of which is that the D~ tag is the daughter of a D*~ and the other that the D* s + decays 
into jDf, with the subsequently decaying into /i + z/. The kinematical constraints, in the 
e + e~ center-of-mass frame, are 

PD a + PD* = (7) 

Eqm = E Ds + E D * 

E CM , M 2 Dt -M 2 Ds E CM M 2 Dt -M 2 Ds 



2 2£'cm 2 2i?cM 

Mu : - M Ds = 143.8 MeV. 

In addition, we constrain the invariant mass of the D~ tag to the known D s mass. This 
gives us a total of 7 constraints. The missing neutrino four-vector needs to be determined, 
so we are left with a three-constraint fit. We perform a standard iterative fit minimizing x 2 ■ 
As we do not want to be subject to systematic uncertainties that depend on understanding 
the absolute scale of the errors, we do not make a x 2 cut but simply choose the photon and 
the decay sequence in each event with the minimum x 2 ■ 

In this analysis, we consider three separate cases: (i) the track deposits < 300 MeV in the 
calorimeter, characteristic of a non-interacting pion or a muon; (ii) the track deposits > 300 
MeV in the calorimeter, characteristic of an interacting pion, and is not consistent with being 
an electron; (iii) the track satisfies our electron selection criteria; for a track to be called an 
electron, we require that the momentum measurement in the tracking system and the energy 
deposited in the Csl calorimeter are close to being equal within errors, and we also require 
that dE/dx and RICH information be consistent with expectations for an electron. Then 
we separately study the MM 2 distributions for these three cases. The separation between 
muons and pions is not complete. Case (i) contains 98.8% of the muons but also 55% of the 
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pions, while case (ii) includes 1.2% of the muons and 45% of the pions [30|. Case (iii) does 
not include any signal but is used later to search for — ► e + v. For cases (i) and (ii) we 
insist that the track not be identified clS cL kaon. 

C. The Expected MM 2 Signal Spectrum 

For the /i + z/ final state the MM 2 distribution can be modeled as the sum of two Gaussians 
centered at zero (see Eq. (j3J). A Monte Carlo simulation of the MM 2 is shown in Fig. [HI A fit 
using the two-Gaussian shape gives o\ = 0.0240 GeV 2 , a 2 = 0.0851 GeV 2 , / = 0.275, which 
results in a = 0.0346 ± 0.0002 GeV 2 . 
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FIG. 8: The MM 2 resolution from Monte Carlo simulation for Df — > The curve is the sum 

of two-Gaussians with means constrained to be the same. 

We check the resolution using data. The mode Df — ► T(°K + provides an excellent testing 
ground. 1 We search for events with at least one additional track identified as a kaon using 
the RICH detector, in addition to a Dj tag. We allow events with no more than two other 
additional charged tracks, to allow for the presence of K° decays, and we do not apply the 
greater than 300 MeV extra energy cut. The MM 2 distribution is shown in Fig. [9j Fitting 
this distribution to a two-Gaussian signal shape gives a MM 2 resolution a = 0.0338 ± 0.0014 
GeV 2 in agreement with Monte Carlo simulation which gives 0.0344±0.0003 GeV 2 . (The 
backgrounds are discussed in section HVl ) We note that the resolution here is larger than in 
our previous work [l[ . This is mainly due to the use of additional decay modes with photons, 
and enlarging of the solid angle for candidate muons. 



1 In this paper the notation K K + refers to the sum of K K + and K°K + final states. 
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FIG. 9: The MM 2 distribution for events with an identified K + track. The kinematic fit has been 
applied. The data are shown as points with error bars. The long-dashed curve shows the calculated 
yield of r/K + events. The solid curve shows the results of a fit to the data, where the dotted curve 
is the sum of two Gaussians centered at the square of the K° mass, and the dashed and dot-dashed 
lines refer to the sideband, and combinatoric backgrounds, respectively. 

For the r + z/, r + — > n + V final state a Monte Carlo simulation of the MM 2 spectrum is 
shown in Fig. [101 The extra missing neutrino results in a smeared distribution. 

D. MM 2 Spectra in Data 

The MM 2 distributions from data are shown in Fig. [TTJ The overall signal region we 
consider is —0.1 < MM 2 < 0.20 GeV 2 . The higher limit is imposed to exclude background 
from 7]7T + and K°tt + final states. There is a clear peak in Fig. flTT a) due to — > fi + u. 
Furthermore, the region between the /z + z/ peak and 0.20 GeV 2 has events that we will show 
are mainly due to the — > r + u, r + — > tt + v decay. The events in Fig. [TTT b) below 0.20 
GeV 2 are also mostly due to t + v decay. 

E. Background Evaluations 

We consider the background arising from two sources: one from real decays and the 
other from the background under the single-tag signal peaks. For the latter, we obtain the 
background from data by using a two-dimensional extended likelihood fit in invariant mass 
and MM 2 . 

The background from real decays is studied by identifying each possible source mode 
by mode. For the /i + z/ final state, the only possible background within the signal region is 
— ► 7r + 7T°. (Recall that any such events are also heavily suppressed by the extra photon 
energy veto of 300 MeV.) We search for this mode by examining the ti + it invariant mass 
spectrum in events where we have selected a D~ tag based on invariant mass only. See 
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Fig. [12j There is no peak at the mass. Fitting to a linear background plus a Crystal- 
Ball signal function, whose shape is fixed by Monte Carlo gives -3.6±8.4 events. Setting the 
mean value to zero results in an upper limit B(Df — > 7r + 7r°) < 3.8 x 10~ 4 at 90% confidence 
level. Multiplying by our 43859 tags and the 2% inefficiency for detecting a > 300 MeV 
photon from the 7r° decays, results in an upper limit of 1/4 of an event in our sample, which 
we ignore. 

For the t + u, t + — > -k + V final state the real backgrounds include, in addition to the 
7T + 7r° background discussed above, semileptonic decays, possible 7r + 7r°7r° decays, other r + 
decays, and small amounts of K°tt + and 7]tt + whose low MM 2 tails leak into the signal re- 
gion. Semileptonic decays involving muons are equal to those involving electrons shown in 
Fig. rmf c) . One event consistent with the electron hypothesis is present. The 7r + 7r°7r° back- 
ground level is estimated by considering the 7r + 7r + 7r~ final state whose measured branching 
fraction is (1.11±0.07±0.04) % [31]. This mode has large contributions from / (980)7r + and 
other 7T + 7r _ resonant structures at higher mass 32l . The 7r + 7r°7r° mode will also have these 



contributions, but the MM 2 opposite to the ir + will be at large mass. The only component 
that can potentially cause background is the non-resonant component measured by FOCUS 
as (17±4)% [32] • This background has been evaluated by Monte Carlo simulation as have 
backgrounds from other r + decays, K°tt + and i]ir + . The backgrounds are enumerated in 
Table HTl We show in Fig. [T3]the sum of all backgrounds and a fit to a quadratic polynomial 
over the MM 2 range of interest. 
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FIG. 11: The MM 2 distributions from data for events with a Dj reconstructed in a tag mode, 
an additional positively charged track and no neutral energy clusters above 300 MeV. (a) Case (i) 
when the single track deposits < 300 MeV of energy in the calorimeter. The peak near zero is from 
Df — > events, (b) Case (ii): the track deposits > 300 MeV in the crystal calorimeter but is 
not consistent with being an electron, (c) Case (iii): the track is identified as an electron. 



TABLE II: Background estimates for the data in the signal region —0.1 <MM 2 < 0.2 GeV 2 . (We 
assume B(D+ -> t + v) = 6.2 ± 0.7%.) 



Final State 


B{%) 


# of events case(i) 


# of events case (ii) 


T + — > 7r + 7T°z7 


1.6±0.2 


2.06±0.34 


1.43±0.36 


T + — > jl + UU 


l.ldbO.l 


1.60±0.24 





D+ -» 7r+7r vr 


1.1 (estimate) 


0.12 


0.12 


Df -» K°TT + 


0.24±0.03 


1.3 ±0.3 


1.1 ±0.3 


D+ -» r/7T + 


1.5 ± 0.2 


1.1 ± 0.3 


0.9 ± 0.3 


Sum 




6.2±0.7 


3.5±0.6 



III. LEPTONIC BRANCHING FRACTIONS 

The result of the two-dimensional unbinned maximum liklihood fit to the sum of the MM 2 
distributions for case (i) and case (ii) is shown in Fig. [141 The other dimension in the fit is 
the invariant mass spectrum. Here we constrain the T + u/fi + u ratio to the Standard Model 
value. From the fit we extract 235.5±13.8 signal events. The efficiency is given by the 
product of the tracking and particle identification efficiencies, equal to 86.7%, the maximum 
extra photon energy cut of 300 MeV, equal to 98.7% and the tag bias of 105.2%, giving an 
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FIG. 12: The invariant tt + tt mass. The curves show results of a fit using a linear background 
(solid) plus Gaussian signal function (dashed), where the width of the Gaussian is fixed to a value 
determined by Monte Carlo simulation. 
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FIG. 13: The background rates for real decays from Monte Carlo as a function of MM 2 for 
case (i) and case (ii). The data are fit to quadratic polynomials. 

overall efficiency of 90.0%. Using our tag sample of 43859±936 events, we find an effective 
branching ratio 

B eS (Dj -> fi+v) = (0.597 ± 0.037 ± 0.017)% (8) 
The radiative correction of 1% reduces this to 

B eS {D+ -> n+v) = (0.591 ± 0.037 ± 0.018)% (9) 

(From now on we will only quote radiatively corrected results in this paper.) This is our 
most accurate result within the context of the Standard Model. 

We can also analyze the data by not constraining the t + v/^i + v ratio. We then fit the 
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FIG. 14: The results of the two-dimensional fit to the sum of case (i) and case (ii) data. The 
data are shown as points with error bars, (a) The projection of the invariant mass distribution; 
the straight dashed line shows the background while the curve is the sum of the background and 
a two-Gaussian signal function, (b) The projection of the MM 2 distribution; the dotted (black) 
curve is the two-Gaussian signal function for [i + v, the long-dashed (purple) curve shows the t + u, 
t + — > 7r + I7 signal, the dashed (red) line shows the background from non-D~ events below the signal 
peak, while the dot-dashed (green) curve shows the background from real Df events. The solid 
(blue) curve represents the sum of all contributions. 



case (i) distribution and find 

B(D+ -> = (0.565 ± 0.045 ± 0.017)% (10) 

The results of the fit are shown in Fig. [15j 

By simultaneously fitting the case (i) and case (ii) distributions, constraining the ratio 
of events to be 98.8:1.2 and the ratio of t + v events to be 55:45, in the two cases, re- 
spectively, we find 125.6±15.7 t+u, r+ -> n + V events. Using B(r + -> tt+z/) = (10.90±0.07)% 
12|. we measure 

B{D+ tV) = (6.42 ± 0.81 ± 0.18)% (11) 



The fit results are shown in Fig. 

The systematic errors on these branching fractions are given in Table IIHI The error on 
track finding is determined from a detailed comparison of the simulation with double tag 
events where one track is ignored. The particle identification on the fi + track arises from 
the fact that we veto kaons. The error on the photon veto efficiency, due to the 300 MeV/c 
extra shower energy cut, is determined using double tag events where both D~ and are 
reconstructed in the same modes that we use for tagging. What we measure is a product of 
the efficiencies for two events each with different tags. This is translated into an efficiency 
for a single tags. The error on the number of tags ±2% is assigned by varying the fitting 
functions and ranges. In addition there is a small error of ±0.6% on the t + u branching 
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FIG. 15: The results of the two-dimensional fit to the case (i) data. The data are shown as points 
with error bars, (a) The projection of the invariant mass distribution; the straight dashed line shows 
the background while the curve is the sum of the background and a two-Gaussian signal function, 
(b) The projection of the MM 2 distribution; the dotted (black) curve is the two-Gaussian signal 
function for n + v, the long-dashed (purple) curve shows the r + ^, r + — ► tt + V signal, the dashed 
(red) line shows the background from non-D~ events below the signal peak, while the dot-dashed 
(green) curve shows the background from real events. The solid (blue) curve represents the 
sum of all contributions. 



fraction due to the uncertainty on the r + decay fraction to n + V. Additional systematic 
errors arising from the background estimates are at the 1% level. The error on the radiative 
correction is taken as 100% of its value of 1%. When we use only one of the two cases to find 
a result an additional 1% error is included due to the minimum ionization discrimination of 
300 MeV in the calorimeter. 



TABLE III: Systematic errors on determination of the Df — > n + v branching fraction. 



Error Source Size (%) 

Track finding 0.7 
Particle identification of fi + 1.0 

MM 2 width 0.2 

Photon veto 0.4 

Background 1.0 

Number of tags 2.0 

Tag bias 1.0 

Radiative Correction 1.0 

Total 31) 
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FIG. 16: The MM 2 distributions from the two-dimensional fit done simultaneously to (a) case (i) 
and (b) case (ii) data. The data are shown as points with error bars. The dotted (black) curve is 
the two-Gaussian signal function for the long-dashed (purple) curve shows the t + i>, t + tt + V 
signal, the dashed (red) line shows the background from non-Z?" events below the signal peak, 
while the dot-dashed (green) curve shows the background from real Df events. The solid (blue) 
curve represents the sum of all contributions. 
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Lepton universality in the Standard Model requires that the ratio R from Eq. (j2j) be equal 
to a value of 9.76. We measure 



rm+ -> t + u) 

fl " fb^7M = 1L4±L7±a2 - (12) 

Here the systematic error is dominated by the uncertainty on the minimum ionization cut 
that we use to separate the /x + z/ and t + v regions at 300 MeV. We take this error as 2%, since 
a change here affects both the numerator and denominator. The ratio is consistent with the 
Standard Model prediction. Current results on D + leptonic decays also show no deviations 

a. 

We also measure the CP violating asymmetry. The fit constraining the r + u/fi + u ratio 
to the SM value of 9.76 yields 124.5±9.9 fj+u and 110.8±9.6 \Tv events. We also find 
21807±581 D~ tags and 21370±581 D+ tags. Then 



r(D+ ->/*+»/) -r(fl7->/z-i/) 



0.048 ±0.061, (13) 



showing no evidence of CP violation. 

The one detected electron opposite to our tags allows us to set an upper limit of 

B(Df -> e + v) < 1.2 x 10" 4 (14) 

at 90% confidence level; this is also consistent with Standard Model predictions and lepton 
universality. 



IV. CHECKS OF THE METHOD 

We perform an overall check of our procedures by measuring B(D^ — > TC°K + ) which 
has been previously determined. For this measurement we compute the MM 2 [Eq. (jSJ)] using 
events with an additional charged track but here identified as a kaon. These track candidates 
have momenta of approximately 1 GeV/c; here our RICH detector has a pion to kaon fake 
rate of 1.1% with a kaon detection efficiency of 88.5% [25]. For this study, we do not veto 
events with neutral energy deposits >300 MeV, or with less than three additional tracks 
beyond the tag, because of the presence of the K°. 

Events from the r]n + mode where the 7r + fakes a K + are very rare and would not peak 
at the proper MM 2 . The mode 7]K + does contribute a background at a somewhat higher 
MM 2 of 0.30 GeV 2 and causes a small asymmetric tail on the high side of the peak. The 
branching fraction measured by CLEO is (0.14±0.03)% [34]. We predict a total of 47±10 
events from this source, that we include as a fixed component in our fit. 

We perform the same two-dimension fit in invariant mass and MM 2 as used for the 
signal, except that here the MM 2 is calculated with respect to the K + hypothesis, and an 
additional 7]K + component is added. The MM 2 distribution for events in the signal MM* 2 
region is shown in Fig. O The peak near 0.25 GeV 2 is due to the decay mode of interest. 
The backgrounds are the same as defined for the fi + u distributions above. 

The fit yields 1036±41 events. In order to compute the branching fraction we use the effi- 
ciency of detecting the kaon track, 77.0%, including radiation [35j, the particle identification 
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efficiency of 88.5%, and take into account that it is easier to detect tags in events containing 
a K°K + decay than in the average D S D* event due to the track and photon multiplicities, 
which gives a 3% correction. 2 These rates are estimated by using Monte Carlo simulation. 
We determine 

B(D+ -> K°K+) = (3.06 ± 0.14 ± 0.09)%, (15) 

where the systematic errors are listed in Table IIV1 We estimate the error from the signal 
shape by taking the change in the number of events when varying the signal width of the 
two-Gaussian function by ±lcr. The error on the background shapes is given by varying the 
shape of the background fit. The error on the particle identification efficiency is measured 
using two-body D° decays [25[ . The other errors are the same as described in Table IIHI 
Again, the largest component of the systematic error arises from the number of tag events 
(2%). In fact, to use this result as a check on our procedures, we need only consider the 
systematic errors that are different here than in the fi + u case. Those are due only to the 
signal and background shapes, the i]K + contribution and the particle identification cut. 
Those systematic errors are small. 

To determine absolute branching fractions of charm mesons, CLEO-c uses a method where 
both particles are fully reconstructed (so called "double tags" ) and the rates are normalized 
using events where only one particle is fully reconstructed. Our result using this method for 
Bm+ -> K° S K+) = (1A9±0.07 ± 0.05)%, which when doubled becomes (2.98±0.14 ± 0.10)% 
[3l|. This is in excellent agreement with the number in Eq. (fl~5l) . These results are not 
independent. 



TABLE IV: Systematic errors on determination of the Df — > K K + branching fraction. 



Error Source 


Size (%) 


Track finding 


0.7 


Particle identification of /x + 


1.0 


r/K + branching fraction 


0.6 


MM 2 width 


0.2 


Background 


1.0 


Number of tags 


2.0 


Tag bias 


1.0 


Total 


2.8 



We also performed the entire analysis on a Monte Carlo sample that corresponds to an 
integrated luminosity 8 times larger than the data sample. The input branching fraction in 
the Monte Carlo is 0.61% for /i + z/ and 5.99% for t + u, while our analysis measured B (D£ — > 
/i + z/) =(0.607±0.013)% for fi + u and t + v combined. The individual rates are B(Df — > 
/i+z/) =(0.615±0.016)%, and B(D+ -> t+v) =(6.02±0.27)%. 



2 The tag bias is less here than in the [i + v case because of the K° decays and interactions in the detector. 
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V. THE DECAY CONSTANT AND CONCLUSIONS 



Using our most precise value for B(Df — > [i^v) from Eq. ([8]), that is derived using both 



our n + v and t + u samples, and Eq. ([I]) with a D s lifetime of (500±7) x 10 15 s 12J, we extract 



f D + = 263.3 ± 8.2 ± 3.9 MeV. (16) 

This result has been radiatively corrected. These results supersede all our previous mea- 
surements of the -D+ — > and — > t + u, t + — > n + u branching fractions which use data 
samples that are subsumed in this paper. Using the CLEO-c result based on an analysis of 
D+ -> t+v, r+ -> e+vi? 0, of f D + = (252.5 ± 11.1 ± 5.2) MeV, we derive a CLEO-c average 
value of 

f D + = 259.5 ± 6.6 ±3.1 MeV. (17) 
We combine with our D + result 0] 

f D + = 205.8 ± 8.5 ± 2.5 MeV (18) 

and find a value for 

^ = 1.26 ±0.06 ±0.02, (19) 

JD+ 

where only a small part of the systematic error cancels in the ratio of our two measurements. 
Our new measurements are compared with other measurements in Table |V] 



TABLE V: Our results for B{Df — > B(Df — » t + v), and f D + compared with previous 

measurements. Results have been updated for the new value of the D s lifetime of 0.5 ps 
ALEPH combines both measurements to derive a value for the decay constant. (This table adopted 
from Table I of ref. (7|.) 



Exp. 


Mode 


B 


(%) 


f D + (MeV) 


CLEO-c 


A 1 




(5.65 ±0.45 ±0.17) • 10~ 3 




257.3 ±10.3 ±3.9 


CLEO-c 






(6.42 ±0.81 ±0.18) • 10~ 2 




278.7 ± 17.1 ±3.8 


CLEO-c 






combined above 2 results using SM 




263.3 ±8.2 ±3.9 


CLEO-c [2] 






(5.30 ±0.47 ±0.22) • 10~ 2 




252.5 ± 11.1 ± 5.2 


CLEO-c 






combined all CLEO-c results 




259.5 ±6.6 ±3.1 


Belle* [21] 


yr 




(6.38 ±0.76 ±0.52) • 10~ 3 




274 ± 16 ± 12 


Average of CLEO and Belle results above, radiatively corrected 




261.2 ±6.9 


CLEO [14] 






(6.2 ±0.8 ±1.3 ±1.6) -10~ 3 


3.6±0.9 


273 ± 19 ± 27 ± 33 


BEATRICE [15] fi~ 




(8.3 ±2.3 ±0.6 ±2.1) • 10~ 3 


3.6±0.9 


312 ± 43 ± 12 ±39 


ALEPH [16] 




"i/ 


(6.8 ± 1.1 ± 1.8) • 10~ 3 


3.6±0.9 


282 ± 19 ± 40 


ALEPH [16] 




"1/ 


(5.8 ±0.8 ±1.8) • 10~ 2 






L3 [17] 




~v 


(7.4 ± 2.8 ± 1.6 ± 1.8) • 10~ 2 




299 ± 57 ± 32 ± 37 


OPAL [18] 


t 1 


~v 


(7.0 ±2.1 ±2.0) • 10" 2 




283 ± 44 ± 41 


BaBar [19] 




-v 


(6.74 ± 0.83 ± 0.26 ± 0.66) • 10~ 3 


4.71±0.46 


283±17±7±14 



f This result has been radiatively corrected by multiplying the measured branching ratio by 99%. 
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Most measurements of — > £ + z/ are normalized with respect to B(Df — > </>7r + ) = B^-. 
An exception is the OPAL measurement which is normalized to the D s fraction in Z° events 
that is derived from an overall fit to heavy flavor data at LEP [33]. It still, however, relies on 
absolute branching fractions that are hidden by this procedure, and the estimated error on 
the normalization is somewhat smaller than that indicated by the error on B^ available at 
the time of their publication. The L3 measurement is normalized taking the fraction of D s 
mesons produced in c quark fragmentation as 0.11±0.02, and the ratio of D*/D s production 
of 0.65±0.10. The ALEPH results use B^ for their /i + z/ results and a similar procedure as 
OPAL for their t + v results. We note that the recent BaBar result uses a larger B^ than 
the other results. The CLEO-c determination of f D + using the modes in this paper is the 
most accurate to date. 

f 4 

Theoretical models that predict f D + and the ratio y^- are listed in Table IVT1 Upper limits 

TABLE VI: Theoretical predictions of f D +, and f D +/ Jd+- QL indicates quenched lattice 

calculations. (This table adopted from Table II of ref. [7|.) 



Model 


f D f (MeV) 


f D+ (MeV) 




Lattice (HPQCD+UKQCD) [6j 


241 ±3 


208 ±4 


1.162 ±0.009 


Lattice (FNAL+MILC+HPQCD) [5] 


249 ± 3 ± 16 


201 ± 3 ± 17 


1.24 ±0.01 ±0.07 


QL (QCDSF) [37] 


220 ± 6 ± 5 ± 11 206 ± 6 ± 3 ± 22 1.07 ± 0.02 ± 0.02 


QL (Taiwan) [38] 


266 ± 10 ± 18 


235 ± 8 ± 14 


1.13 ±0.03 ±0.05 


QL (UKQCD) [39] 


236 ± 8±\l 


210 ± 10±\l 


1.13 ±0.02^02 


QL [40] 


231 ± 1215 


211 ± U±l 2 


1.10 ±0.02 


QCD Sum Rules [41] 


205 ± 22 


177 ±21 


1.16 ±0.01 ±0.03 


QCD Sum Rules [42] 


235 ± 24 


203 ± 20 


1.15 ±0.04 


Field Correlators [43] 


210 ±10 


260 ± 10 


1.24 ±0.03 


Quark Model [44] 


268 


234 


1.15 


Quark Model [45] 


248±27 


230±25 


1.08±0.01 


LFQM (Linear) [46J 


211 


248 


1.18 


LFQM (HO) [46] 


194 


233 


1.20 


LF-QCD [47] 


253 


241 


1.05 


Potential Model [48] 


241 


238 


1.01 


Isospin Splittings [49] 




262 ± 29 





on f D + and fjj a of 230 and 270 MeV, respectively, have been determined using two-point 
correlation functions by Khodjamirian [36]. Our result for fu a is higher than most theoretical 
expectations. The average of this new f D + result with the CLEO-c result based on the decay 
mode D+ -> t + v, r+ -> e + W, is (259.5 ± 7.3) MeV (see Table [VI). This rate differs by 2.3 
standard deviations from the Follana et al. prediction [6], and is close to, but does not 
saturate the Khodjamirian bound. (Averaging in the Belle result for — ► which 
is also an absolute branching fraction measurement raises the difference to 2.6 standard 
deviations.) If the difference with Follana et al. were to persist with the advent of more 
precise measurements it could be explained by a deficiency in the Lattice calculation, or 
physics beyond the Standard Model. (We note that both unquenched LQCD calculations 
agree with CLEO's result for fo+-) Either of the two leptoquark models of Dobrescu and 
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Kronfeld could explain a discrepancy in the Dj~ c ase [8J. They also have a charged Higgs 
model. Hewett and Akeroyd and Chen ll[ pointed out that leptonic decay widths 
are modified by new physics. Specifically, for the D + and Df, in the case of the two-Higgs 
doublet model (2HDM), Eq. ([!]) is modified by a factor r q multiplying the right-hand side: 



(20) 



where is the charged Higgs mass, Mjj is the mass of the D meson (containing the light 
quark q) , m c is the charm quark mass, m q is the light-quark mass, and tan (3 is the ratio of 
the vacuum expectation values of the two Higgs doublets. (Here we modified the original 
formula of 11] to take into account the charm quark coupling [loj].) To get an enhancement 
in the rate, the m c term must be inducing the effect, which implies that both the D + and 
would see a similar effect in contradiction to the trends in our data. Another explanation 
in a model based on R-parity violating supersymmetry has been given by Kundu and Nandi 
51 . 
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